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Abstract 

We discover multi-Hamiltonian structure of complex Monge- Ampere 
equation (CM A) set in a real first-order two-component form. There- 
fore, by Magri's theorem this is a completely integrable system in four 
real dimensions. We start with Lagrangian and Hamiltonian densities 
and obtain a symplectic form and the Hamiltonian operator that de- 
termines the Dirac bracket. We have calculated all point symmetries 
of two-component CM A system and Hamiltonians of the symmetry 
flows. We have found two new real recursion operators for symme- 
tries which commute with the operator of a symmetry condition on 
solutions of the CM A system. These operators form two Lax pairs 
for the two-component system. The recursion operators, applied to 
the first Hamiltonian operator, generate infinitely many real Hamil- 
tonian structures. We show how to construct an infinite hierarchy of 
higher commuting flows together with the corresponding infinite chain 
of their Hamiltonians. 

PACS numbers: 04.20.Jb, 02.40.Ky 

AMS classification scheme numbers: 35Q75, 83C15 

1 Introduction 

In earlier paper [1] we presented complex multi-Hamiltonian structure of 
Plebanski's second heavenly equation [2], which by Magri's theorem [3] proves 
that it is a completely integrable system in four complex dimensions. We ex- 
pect that Plebanski's first heavenly equation also admits multi-Hamiltonian 
structure, since these two equations, governing Ricci-flat metrics with (anti- 
) self-dual Riemann curvature 2- form, are related by Legendre transformation 



of corresponding heavenly tetrads [2]. However, since both Plebanski's equa- 
tions are complex, their solutions are potentials of the complex metrics that 
satisfy Einstein equations in complex four-dimensional spaces. In the case of 
complex Monge- Ampere equation (CM A), that governs (anti-) self-dual grav- 
ity in real four-dimensional spaces with either Euclidean or ultra-hyperbolic 
signature, we have an additional condition that symplectic, Hamiltonian, and 
recursion operators all should be real. Furthermore, the transformation be- 
tween the two heavenly equations cannot be applied to transform the second 
heavenly equation to CM A because the latter equation is real and therefore 
the multi-Hamiltonian structure of CM A cannot be obtained by transform- 
ing the multi-Hamiltonian structure of the second heavenly equation given 
in [1]. Therefore, in this paper we analyze the complex Monge- Ampere equa- 
tion independently of our previous work and obtain real recursion operators 
for symmetries and real multi-Hamiltonian structures of CM A. 

In section [2] we start with the complex Monge-Ampere equation in a 
two-component first-order evolutionary form with the Lagrangian that is 
appropriate for Hamiltonian formulation. In section [3] we present a sym- 
plectic structure and Hamiltonian structure of this CM A system. In section 
H] we transform the Hamiltonian density and Hamiltonian operator to real 
variables and introduce a convenient notation needed later to arrive at a 
compact form of recursion operators and higher Hamiltonian operators. In 
section we derive a symmetry condition, that determines symmetries of 
CM A system, in a two-component form and real variables, using our new 
notation. We have calculated all point symmetries of the CM A system and 
Hamiltonians of the symmetry flows that yield conservation laws for the 
CM A system. In section [6] we obtain two new real recursion operators for 
symmetries which commute with an operator of the symmetry condition on 
solutions of the CM A system. Moreover, these two couples of operators form 
two Lax pairs for the two-component system. The recursion operators, being 
applied to the first Hamiltonian operator, yield further Hamiltonian struc- 
tures and bi-Hamiltonian reresentations of the CM A system. Repeating this 
procedure for the second Hamiltonian operator, we could generate infinitely 
many Hamiltonian structures of the CM A system. This multi-Hamiltonian 
structure of the CM A system proves its complete integr ability in the sense 
of Magri and hence complete integrability of the (anti-) self-dual gravity in 
four real dimensions with either Euclidean or ultra-hyperbolic signature. In 
section [7] we construct an infinite hierarchy of higher flows and show the way 
of calculating a corresponding infinite chain of higher Hamiltonians. 
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2 Complex Monge- Ampere equation in first- 
order evolutionary form and its Lagrangian 

Four- dimensional hyper-Kahler metrics 

ds 2 = undz dz + u^dz^dz + U2idz 2 dz x + u 2 2~dz 2 dz 2 (2.1) 

satisfy Einstein field equations with either Euclidean or ultra-hyperbolic 

signature, if the Kahler potential u satisfies elliptic or hyperbolic complex 
Monge- Ampere equation 

M11M22 - = £ (2.2) 

with e = ±1 respectively [2]. Here u is a real- valued function of the two 
complex variables z x ,z 2 and their conjugates z 1 ,^ 2 , the subscripts denoting 
partial derivatives with respect to these variables. Such metrics are Ricci-flat 
and have (anti-)self-dual curvature. 

In order to discuss the Hamiltonian structure of CM A (12.21) . we shall 
replace the complex conjugate pair of variables z 1 , z 1 by the real time variable 
t = 2^R.z 1 and the real space variable x = 2QZ 1 and change the notation for 
the second complex variable z 2 = w. Then (12. 2p becomes 

\^Ht ^xx^^ww ^tw^tui 'U'xw'U'xw ^i^tw^xw ^xw^tui) (2-3) 

Now we can express ( 12.31) as a pair of first-order nonlinear evolution equations 
by introducing an auxiliary dependent variable v = u t 

{u t = v 
1 / _ x (2.4) 

so that finally (12.21) is set in a two-component form. For the sake of brevity 
we shall henceforth refer to ( 12.41) as CM A system. 

The Lagrangian density for the original form (I2.2p of the complex Monge- 
Ampere equation was suggested in [4] 

L = - \u x u\u 2 2 + U2U2UH - U1U2U21 - U2U1U12} + eu, (2.5) 
o 

but this must be cast into a form suitable for passing to a Hamiltonian. This 
requires that the form of a Lagrangian should be appropriate for applying 
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Dirac's theory of constraints [5]. We choose the Lagrangian density for the 
first-order CM A system (12 .4p to be degenerate, that is, linear in the time 
derivative of unknown Ut and with no Vt- 

6 

+ u t (2i(u w u xw - UyjU xw ) + 6toot))} + eu (2.6) 

which, after substituting v = u t , coincides with our original Lagrangian (12.51) 
up to a total divergence. 



3 Symplectic and Hamiltonian structures 

Since the Lagrangian density (12.61) is linear in u t and has no v t , the canonical 
momenta 

dL i 

T^u ~7{ TlK'U'w'U'X'w V'w'U'xw) V^ ww 

out 3 

7T V = ^ = (3.1) 

dv t 

cannot be inverted for the velocities u t and v t and so the Lagrangian is 
degenerate. Therefore, according to the Dirac's theory [5], we impose them 
as constraints 

(fiu 7T« T^yU'W^xw ^w^xw) VU^jyi 

<f> v = tc v = (3.2) 

and calculate the Poisson brackets of the constraints (more details of the 
procedure were given in [1]) 

K ik = [<f>i(x, w, w), (f>k(x', w', w)] (3.3) 

collecting results in a 2 x 2 matrix form, where the subscripts run from 1 
to 2 with 1 and 2 corresponding to u and v respectively. This yields the 
symplectic operator K that is the inverse of the Hamiltonian operator Jo: 

| ipw i^xw)D w -\- {y w -\- iu X w)D^[, -\- v W yj Uww j ^ ^ 
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as an explicitly skew-symmetric local operator. A symplectic 2-form is a 

volume integral ft = J udxdwdw of the density 

v 

1-1 1 
uj = - du l AKij du 3 = -(vu,—iu xii j)duAdu w +-(v w +iu xw )duAdu i jj+u W u ! dvAdu 

2 2 2 (3-5) 

where u l = u and u 2 = v. In u, under the sign of the volume integral, we 
can neglect all the terms that are either total derivatives or total divergencies 
due to suitable boundary conditions on the boundary surface of the volume. 
For the exterior differential of this 2-form we obtain 

du = —idu x A du w A duu, = — (i/3)(^D x (du A du w A duw) (3.6) 

+ D w {du x A du A du^) + D^[du x A du w A du)^j 

that is, a total divergence which is equivalent to zero, so that the 2-form Q 
is closed and hence symplectic. The Hamiltonian operator J is obtained by 
inverting K in (13.41) 

Jo = (3.7) 

/ -J- 

\ 2m 2 w w ^h? + 2t? * ™ 2?? 

that is explicitly skew- symmetric. It satisfies the Jacobi identity due to (13. 6p . 
The Hamiltonian density is 



H x = ir u u t + ir v v t - L 



with the result 



Hi = - (3v - U x ^jU ww r U> w U- w 'U xx ~\~ 'U X (^U- W / U XW ~\~ U-w'U'xw^l Sll. (3.8) 
L J 

CM A system can now be written in the Hamiltonian form with the Hamil- 
tonian density Hi defined by (13. 8|) 

(:H(S) ™ 

where 5 U and 6 V are Euler-Lagrange operators [6] with respect to u and 
v applied to the Hamiltonian density Hi (they correspond to variational 
derivatives of the Hamiltonian functional / HidV). 
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4 Transformation to real variables 



In the case of CMA, that governs (anti-)self-dual gravity with either Eu- 
clidean or ultra-hyperbolic signature, we have an additional condition that 
all the objects in the theory, in particular a recursion operator, should be 
real. Therefore, we transform the Hamiltonian density together with the 
symplectic and Hamiltonian operators to the real variables y = 2^tw and 
z = 29m The Hamiltonian density in the real variables becomes 
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(3v 2 - u 2 x )A(u) - (u 2 + u 2 z )u xx + 2u x (u y u xy + u z u xz ) 



eu. 



where A{u) = u yy + u zz , which simplifies after cancelling terms that are total 
derivatives to 



Hi =2 b 2A ( M ) - u xx{u 2 y + U 2 Z )} - EU. 



The transformation of the Hamiltonian operator Jq in (13.71) yields 



/ 



Jo 







1 

a 



\ 



V a a 
where we introduce the notation 



1 ^(cDy - bD z ) + (D y c - D z b)± 



a = A(u), b 



V'xy V Z1 C Vy -\- U XZ} 



Q 



b 2 + c 2 + e 



(4.1) 



(4.2) 



(4.3) 



that we will use from now on throughout the paper, with D y , D z designating 
operators of total derivatives with respect to y, z respectively and A = D 2 + 
D 2 is the two-dimensional Laplace operator. 

The symplectic operator (13.41) in the real variables becomes 



K = J, 1 



cD y - bD z + D y c - D z b 



-a 




(4.4) 



in an explicitly skew-symmetric form. 

CMA system (12.41) in the real variables becomes 



i't 



Jo 



SuHx 
S v Hx 



v 

Q Hxx 



(4.5) 



or u t = v, v t = Q - u xa 
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The four-dimensional hyper-Kahler metrics fl2.ll) in the real variables in 
the notation (14. 3p become 



ds 



Q(dt 2 + dx 2 



+ a(dy 2 + dz 2 ) 



c(dtdy + dxdz) — b(dtdz 



dxdy) 
(4.6) 

The metrics (14. 6 p satisfy Einstein field equations with either Euclidean or 
ultra-hyperbolic signature, if the two- component potential (u, v) in the def- 
initions (14.31) of a, b, c, and Q satisfies the Hamiltonian CM A system (14.51) 
with e = +1 or £ = - 1 respectively . These metrics are again Ricci-flat and 
have (anti-) self-dual curvature. 



5 Symmetries and integrals of motion 

Now, consider Lie group of transformations of the system (14. 5 p in the evolu- 
tionary form, when only dependent variables are transformed, and let r be 
the group parameter. Then Lie equations read 

u T = ip, v T — if) (5.1) 

where $ = ^ ^ ^ is a two-component symmetry characteristic of the system 

(14.51) . The differential compatibility conditions of equations (14. 5 p and (15.11) 
in the form u tT — u Tt = and v tT — v Tt = result in the linear matrix equation 

A($) = (5.2) 

where A is the Frechet derivative of the flow ( 14.51) 

/A -1 \ 

A= \Dl- l{cD z + bD y )D x + § A, D t - l(cD y - bD z ) ) (5 ' 3) 

where the first row of (I5.2p yields (p t — ip. 

Using the software packages LIEPDE and CRACK by T. Wolf [7], run 
under REDUCE 3.8, we have calculated all point symmetries of CM A system 
(I4.5I) . a class of solutions of the matrix equation (I5.2p . We list their generators 
and two-component symmetry characteristics [6], the latter denoted by ip u , (p v 

X x = td t + xd x + ud u , ipl = u-tv- xu x , (p\ = t(u xx - Q) - xv x 
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X 2 = zd y - yd z , ip 2 = yu z - zu y , tp\ = yv z - zv y 
X 3 = d z , (fl = u z , <pl = v z 

X 4 = d y , ifl = u y , (fil = v y (5.4) 
^5 = yd y + zd z + ud u + vd v , Lfl = u-yu y - zu z , (p v 5 = v - yv y - zv z 
X a = a(t,x,y,z)d u + a t {t,x,y,z)d v , tp u a = a, <p v a = a t 
Xp = p z (y, z)d x - P v (y, z)d t , <p% = f3 y v - p z u x , (p° p = (3 y {Q - u xx ) - f3 z v x 

where a(t, x, y, z) is an arbitrary smooth solution of the equations 

A(a) = 0, a tt + a xx = 0, a tz - a xy = 0, a ty + a xz = (5.5) 

whereas /3(y, z) satisfies the two-dimensional Laplace equation A(/3) = 0. 

We shall find the integrals of motion generating the point symmetries 
that serve as Hamiltonians of the symmetry flows 

(:H:H(t*) <™> 

where the symmetry group parameter r plays the role of time for the sym- 
metry flow (15.61) and 7i = Hdxdydz is an integral of the motion along 
the flow f)4.5p . with the conserved density H, that generates the symmetry 
with the two-component characteristic y? u , yv The second equality in (15.61) 
is the Hamiltonian form of Noether's theorem that gives a relation between 
symmetries and integrals. 

We choose here the Poisson structure determined by our first Hamiltonian 
operator J since we know its inverse K given by (14. 4p which is used in the 
inverse Noether theorem 

(&)-*(£) (5 - 7) 

determining conserved densities H corresponding to known symmetry char- 
acteristics (p u ,(p v . 

Using (15.71) . we reconstruct the Hamiltonians of the flows (15.61) for all 
variational point symmetries in (|5.4|) . For the scaling symmetries generated 
by Xi and X$, Hamiltonians do not exist and so they are not variational 
symmetries. For the rotational symmetry generated by X 2 , the Hamiltonian 
is 

H 2 = v{yu z - zu y )A(u) - u x [2(zu y + yu z )u yz + u 2 y + u% (5.8) 
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For the translational symmetries generated by X 3 and X 4 , the corresponding 
Hamiltonians H 3 and H 4 are 



2 
3 



H 3 = vu z A{u) + - u x (u y u zz - u z u yz ) 



2 

Hi = VU V A(U) + - U X (UyUy Z - U Z Uyy). (5.9) 

For the infinite Lie pseudogroups generated by X a and Xp , the Hamiltonians 
are 



1 
2 

and 



H a = avA(u) + - a t (u 2 y + uf) + a(u z u xy - u y u xz ) (5.10) 



Hp = (^-v 2 - (3 z u x vJ A(u) -^u xx (ul + u 2 z ) + ^ul(f3yyUy + (3 yz u z )-s(3yU. 

(5.11) 

In particular, the Hamiltonian of time translations Xp =y = —d t , that is 
Hp =y = Hi, coincides with the Hamiltonian (14.11) of CM A flow. For trans- 
lations in x, Xp =z = d x , the Hamiltonian is Hp =z = —u x vA(u). For a simple 
example of the symmetry X a , X a=z = zd u , the Hamiltonian is H a=z = 
zvA(u) + u x u y , which coincides with the Hamiltonian H (\6.7\i from an in- 
finite chain of Hamiltonians for a hierarchy of higher commuting flows in 
section El 

All of these Hamiltonians of the symmetry flows are conserved densities 
of the CM A flow (1431) . 



6 Recursion operators and bi-Hamiltonian 
representations of CMA system 

Complex recursion operators for symmetries of the heavenly equations of 
Plebanski were introduced in the papers of Dunajski and Mason [8,9]. We 
have used them in our method of partner symmetries for obtaining non- 
invariant solutions of complex Monge-Ampere equation [10, 11] and second 
heavenly equation of Plebanski [11] and, in a two-component form, for gen- 
erating multi-Hamiltonian structure of Plebanski's second heavenly equation 
in [1]. However, for CMA we have an additional condition that the equation 
and its symmetries are real and hence recursion operators should also be real. 
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This condition leads to a couple of real recursion operators in a 2 x 2 matrix 
form. The first one is 



Ri 





QD Z - cD a 



(6.1) 



/ 



V 



D, 



-aD x + bD y + cD. 



DJcD„ - bD. 



D y (cD y 



bD : 



D z (aD x 



bD y - cD, 



-D y a 



-D x D y a 



where A 1 means operator multiplication, and the second recursion operator 
reads 



+ A 



-i 





bD x - QD y c 

( D y {bD z - cD y 



(6.2) 



+ DJ-aD x + bD v + cD z 



D, 



D y (-aD x + bDy + cD z ) + D z (cD y - bD 2 



D y a 
-D, r D y a 



Straightforward, though cumbersome, calculations show that the opera- 
tors R\ and _R 2 commute with the operator A (15. 3p of the symmetry condition 
(15.21) on solutions of equations (14. 5 p and therefore they are indeed recursion 
operators for symmetries of the CM A system. This means that if (<p, ip) is 
obtained by transforming a two-component symmetry characteristic (<p, ip) 
of the system ( 14.51) by the operator R\ or R2 



(6.3) 





where i — 1,2, then (<p,ip) is also a symmetry characteristic of (14.51) . 

Moreover, vanishing of the commutators [i?j,^4], computed without using 
the equations of motion (14. 5p . reproduces the CM A system ( 14. 5p and hence 
the operators Ri and A form two real Lax pairs for the two-component sys- 
tem. Indeed, introducing a short-hand notation F = u t — v, G = v t + u xx — Q, 
$ = F xy — G z and x = Fxz + G y , we rewrite the CM A system in the form 
F — 0, G = 0, so that $ = and x = on its solutions, and the first 
commutator reads 



[Ri,A] 











±[QA(F)-2(M> + cx)]D z -$ 



A 



-1 



x 



(6.4) 
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*{Dj - Dy) - 2 X D y D z 

+ D y A(F)D x - A(F x )D y - A(G)D 

2(D zX D x + D x $D y )D z - D Z A(F)D 



+ [2 Xz - A(G)]D x D y 



D Z A{F) 
D x D y A(F) 



\ 



+ A(G y )D x 



-A(G 



X ) Dy 



and we have a similar expression for [i2 2) >l]. Thus, = implies 

F = 0, G = 0, that is, the CM A system f TCol) . 

These real Lax pairs are formed by the recursion operators for symmetries 
and operator A of the symmetry condition and so they are Lax pairs of the 
Olver-Ibragimov-Shabat type [14,15], which is different from the complex Lax 
pairs suggested by Mason and Newman [12,13] and Dunajski and Mason [8,9] 
and those that we used in [10,11] in relation to partner symmetries, even if we 
set our new Lax pairs in one-component forms. Furthermore, the commutator 
of the complex recursion operator of Mason-Dunajski with the operator of 
the symmetry condition, in a one-component form, reproduces the symmetry 
condition and not the original equation CM A [10]. 

By the theorem of Magri, given a Hamiltonian operator J and a recursion 
operator R, RJ is also a Hamiltonian operator [3]. Thereby, acting by the 
recursion operator R\ on the first Hamiltonian operator Jo (|3.7|) . we obtain 
the second Hamiltonian operator 




D X Dy 



~D X Dy 

DlD z 



+ 



(6.5) 



b 

a 

that is explicitly skew-symmetric. Here Q- = (c 2 — b 2 + e)/a. The proof 
of the Jacobi identity for J\ is lengthy but can be somewhat facilitated by 
using Olver's criterion in terms of functional multi- vectors [6]. 

Similarly, acting by the recursion operator R-i on the Hamiltonian op- 
erator Jo, we obtain a Hamiltonian operator that is another companion for 
Jo- 



J 1 = R 2 J = A 



Dy 

-D^D, 



D X D 
DID 



(6.6) 
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r ^(cD z - aD x ) + (D z c - D x a)-^ + %^D y + D y ^ ) 



\ a 



that is also explicitly skew-symmetric and Q~ = (b 2 — c 2 + e)/a. The Jacobi 
identity for J 1 was also proved by using Olver's criterion in terms of functional 
multi- vectors [6]. 

The flow (14.51) can be generated by the Hamiltonian operator J\ from the 
Hamiltonian density 

H = zvA(u) + u x u y (6.7) 

so that CMA in the two-component form (14.51) admits two Hamiltonian rep- 
resentations 

(:;H(SH(S) (6 - 8 » 

and thus this is a bi- Hamiltonian system. 

The same flow (14.51) can also be generated by the Hamiltonian operator 
J 1 from the Hamiltonian density 

H° = yvA(u) - u x u z (6.9) 

which yields another bi-Hamiltonian representation of the CMA system (14.51) 





= j ' ( 1:"° ) ■ < 6 - io > 

Repeating this procedure n times, we obtain a multi- Hamiltonian repre- 
sentation of the CMA system with the Hamiltonian operators J n = i?" Jo, 
J n = R%J , J"" m = R?R%- m Jo (m = 1, 2, ... n - 1) and corresponding 
Hamiltonian densities. This procedure will be considered in more detail in 
the next section for the operator Ri. Multi-Hamiltonian structure of the 
CMA system proves its complete integrability in the sense of Magri and 
hence the complete integrability of the (anti-) self-dual gravity in four real 
dimensions with either Euclidean or ultra-hyperbolic signature. 

A totally different recursion operator for the (anti-)self-dual gravity in 
complex Einstein spaces was obtained much earlier by Strachan [16] by using 
a Legendre transformed version of the first heavenly equation, that was de- 
rived by Grant [17]. This recursion operator can be factorized which suggests 
a bi-Hamiltonian structure of the resulting evolutionary equation, though 
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that was not completely proved. However, the evolutionary equation and 
the related Hamiltonian structures are expressed in complex variables, with 
the complex "time" t in particular, and with a complex unknown. Therefore, 
the corresponding metric will not correspond to anti-self-dual gravity in real 
Einstein spaces with the Euclidean signature (+ + ++). Furthermore, the 
Poisson bracket contains the unusual operator d^ 1 that could be avoided in 
a two-component formulation. 



7 Infinite hierarchy of higher flows 

The operators Jo and J\ are compatible Hamiltonian operators, i.e. they 
form a Poisson pencil. This means that every linear combination CqJq + C\ J\ 
with constant coefficients C and C\ satisfies the Jacobi identity. This can 
be more easily verified by using the Olver's criterion in terms of functional 
multi-vectors though the calculation is still very lengthy. We know from 
the work of Fuchssteiner and Fokas [18] (see also the survey [19] and ref- 
erences therein) that if a recursion operator has a factorized form, as in 
our case R\ = JiJq 1 = JiK, and the factors Jo and Ji are compatible 
Hamiltonian operators, then R\ is hereditary (Nijenhuis) recursion operator, 
i.e. it generates an Abelian symmetry algebra out of commuting symmetry 
generators. Moreover, Hermitian conjugate hereditary recursion operator 
R\ = Jq X J\ = KJi, acting on the vector of variational derivatives of an 
integral of the flow, yields a vector of variational derivatives of some other 
integral of this flow. Then (16. 8p implies that R\ generates the Hamiltonian 
density H\ from Hq\ 




D x c-D Z Q 
b 



J 1 J 1 



where R\ is defined by 

R\ = 

( (-D x a + D y b + D z c)D y 



+ 



V 



+ (D y c-D z b)D z 
aD z 





(D Z b-D yC )Dy + 

(-D x a + D y b + D z c)D z 
aD x D y 



(7.1) 



(7.2) 



A- 1 . 
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The first higher flow of the hierarchy is generated by J\ acting on the 
vector of variational derivatives of Hi 



where t\ is the time variable of the higher flow. This flow is nonlocal and the 
right-hand side of (17. 3p is too lengthy to be presented here explicitly. 

Now we could generate the next Hamiltonian H 2 of the hierarchy of com- 
muting flows by applying R\ to the vector of variational derivatives of H\\ 

* ( Th\ ) = ( tnl ) = ( fJl ) • ™ 

Therefore, the second higher flow in the hierarchy has a bi-Hamiltonian rep- 
resentation 



where the third Hamiltonian operator J2 is generated by acting with R\ on 
J\\ JiR\ = J\KJ\ = R\J\ = Ji. Acting by J2 on the variational derivatives 
of H , we obtain the relations 

h { S V H ) ~ JlR ' { S V H )- Jl { 5 v Hi ) 

= ( ££ tt ) C") 

where we have used that J\ = J (KJi) = JqR\. From (17. 6p we obtain 
three-Hamiltonian representation of the first higher flow 

(:;:H(t*:H(S)^t£)- <"> 

We could also construct the Hamiltonian H_i such that H is generated 
from H_i by R\ 
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that implies a bi-Hamiltonian representation for the zeroth flow 

(::H(SHK:;)- <™> 

Further we obtain 

*(w:l)= J ^ (*&)=*(*£::) <™> 

and the bi-Hamiltonian representation (16. 8p of the original two-component 
CMA flow becomes a three-Hamiltonian representation of this flow 

We could still continue by applying i?J to the vector of variational deriva- 
tives of H 2 to generate the next Hamiltonian H% 

and obtain a bi-Hamiltonian representation for the next higher flow 
/ u t , \ T I 5 U H 3 \ I 5 U H 2 \ T I 5 U H 2 



where we have used (17.121) and the relation J1KJ1 = R\Ji = J 2 , and so on. 



8 Conclusion 

Our starting point was the symplectic and Hamiltonian structure of the com- 
plex Monge- Ampere equation, set into a two-component evolutionary form. 
We have calculated all point symmetries of the CM A system and also, us- 
ing inverse Noether theorem, Hamiltonians of the flows for all variational 
symmetries. These Hamiltonians yield conservation laws for the CM A flow. 
We have found two real 2x2 matrix recursion operators Ri and R 2 for 
symmetries that commute with the operator A of the symmetry condition 
and hence map any symmetry of the CM A system again into a symme- 
try. The operators R\ and R 2 together with A form two Lax pairs for the 
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two- component CM A system. Acting on the first Hamiltonian operator by 
each recursion operator, we obtain two new Hamiltonian operators according 
to Magri's theorem [3] and two bi-Hamiltonian representations of the com- 
plex Monge-Ampere equation in the two-component form. Repeating this 
action, we could generate an infinite number of Hamiltonian operators and 
hence construct a multi-Hamiltonian representation of the CM A system. We 
show how to construct an infinite hierarchy of higher commuting flows to- 
gether with the corresponding infinite chain of their Hamiltonians by using 
a Hermitian conjugate recursion operator. In particular, we arrive at three- 
Hamiltonian representations for both CM A flow and the first higher flow and 
bi-Hamiltonian representations for the zeroth flow and second higher flow. 
The results of this paper prove complete integr ability of the (anti-) self-dual 
gravity in four real dimensions in the sense of Magri (a multi-Hamiltonian 
representation). 
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